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We present a measurement of the mass difference between t and t quarks in lepton+jets final 
states of ti events in 1 fb _1 of data collected with the DO detector from Fermilab Tevatron Collider 
pp collisions at yfs = 1.96 TeV. The measured mass difference of 3.8 ± 3.7 GeV is consistent with 
the equality of t and t masses. This is the first direct measurement of a mass difference between a 
quark and its antiquark partner. 
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The CPT theorem [l[ , which is fundamental to any lo- 
cal Lorentz-invariant quantum field theory, requires that 
the mass of a particle and that of its antiparticle be iden- 
tical. Tests of CPT invariance for many of the elemen- 
tary particles accommodated within the standard model 



(SM) are available in the literature Despite the fact 
that no violations have ever been observed, it is impor- 
tant to search for the possibility of CPT violation in 
all sectors of the standard model. Because quarks carry 
color, they cannot be observed directly, but must first 
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evolve through quantum chromodynamic (QCD) inter- 
actions into jets of colorless particles. These jet rem- 
nants reflect the characteristics of the initially produced 
quarks, such as their charges, spin states, and masses. 
If the lifetimes of quarks are much longer than the time 
scale for QCD processes, the quarks form hadrons be- 
fore they emerge from collisions, and decay from within 
bound hadronic states. This makes it difficult to measure 
a q — q mass difference because of the model dependence 
of QCD binding and evolution processes. However, since 
the lifetime of the top quark is far shorter than the time 
scale for QCD interactions, the top-quark sector provides 
a way to measure the mass difference less ambiguously Q . 

In this Letter, we report a measurement of the dif- 
ference between the mass of the top quark (t) and that 
of its antiparticle (i) produced in pp collisions at y/s = 
1.96 TeV. Our measurement is based on data correspond- 
ing to ~1 fb _1 of integrated luminosity collected with 
the DO detector 0| during Run II of the Fermilab Teva- 
tron Collider. The events used in this analysis, identical 
to those in Ref. fj|, are top quark pair (tt) events in 
the lepton + jets channel (^+jets) where each top quark 
is assumed to always decay into a W boson and a b 
quark. One of the W bosons decays via W — > Iv into 
two leptons, and the other one through W — > qq' into 
two quarks, and all four quarks (qq'bb) evolve into jets. 

We select events having one isolated electron (muon) 
with transverse momentum pt > 20 GeV and \rj\ < 1.1 
(\rj\ < 2), missing transverse momentum f T > 20 GeV, 
and exactly four jets with pr > 20 GeV and \r)\ < 2.5, 
where the pseudorapidity rj = — In [tan(#/2)], and 6 is 
the polar angle with respect to the proton beam direc- 
tion. At least one of the jets is required to be identi- 
fied as a 6-jet candidate. A minimum azimuthal sepa- 
ration is required between lepton px and j> T vectors to 
further reduce multijet background arising from lepton or 
jet energy mismeasurements. The positively (negatively) 
charged leptons are used to tag the t (?) in each event. To 
reduce instrumental effects that can cause charge depen- 
dent asymmetries in lepton energy scale and resolution, 
solenoid and toroid magnetic field polarities are routinely 
reversed. 

The selected data sample consists of 110 e+jets and 
110 ^+jets events. The W + (W~) boson decays into 
hadrons in 105 (115) events and into leptons in 115 (105) 
events, consistent with invariance under charge conjuga- 
tion. The fraction of it events in this sample is estimated 
to be 74%. The background consists of W^+jets and mul- 
tijet events, with the latter comprising 12% of the entire 
background. 

This analysis uses the matrix element (ME) method 
which relies on the extraction of the properties of the 
top quark (e.g., the mass) through a likelihood technique 
based on probability densities (PD) for each event, cal- 
culated from the ME for the two major processes {it 
and VF+jets production) that contribute to the selected 



£+jets sample. In calculating the PD for it production, 
we include only the leading order (LO) ME from qq — ► it 
production jfjj]. We assume SM-like it production and 
decay, where identical particle and antiparticle masses 
are assumed for b quarks and W bosons but not for top 
quarks. For l^+jets production, we use the ME provided 
in VECBOS 0- The PD for each event is given in terms 
of the fraction of signal (/) and of background (1 — /) in 
the data and the masses of the t (mt) and the t (mj): 

P ovt = A(x)[fP sis (x; m tl m T ) + (1 - f)P hks (x)} , (1) 

where x denotes the measured jet and lepton energies 
and angles, A(x) is a function only of x and accounts 
for the geometrical acceptance and efficiencies, and P s i g 
and Pbkg represent the PD for it and VF+jets production, 
respectively Multijet events are also represented by Pbkg 
since Pbkg 3> Psig for such events Q. 

The free parameters in Eq. [T] are determined from a 
likelihood L(x; m t ,mj, /) constructed from the product 
of the P C vt for all events. Jet energies are scaled by an 
overall jet energy scale (JES) calibration factor derived 
by constraining the reconstructed mass of the two jets 
from W -» qq' decays in it events to 80.4 GeV H %. 
The likelihood is maximized as a function of / for each 
(m,t,mj) hypothesis to determine / best . An integration of 
the likelihood for / = / bcst over the sum m sum = (nit + 
mj)/2 results in a one-dimensional likelihood L(x; A) as 
a function of mass difference A = m t — mj. This is used 
to extract the mean value of A and its uncertainty. A 
similar procedure involving an integration over A gives 
L(x]m sum ) which is used to extract the mean value of 
m sum and its uncertainty. 

The variables in any ME refer to nascent produced par- 
ticles (leptons and partons), but the measured quantities 
correspond to physical leptons and jets. This difference 
is taken into account in the calculation of the event prob- 
ability by convoluting over phase space a transfer func- 
tion, W(y,x), that provides the resolution for the lepton 
in question or a mapping of the observed jet variables in 
an event (x) to their progenitor parton variables (jj): 

p. = _J_ 

Slg ~~ it 
O 

norm 

x J ^da(y;mt,mj)dqxdq2F(q 1 )F(q 2 )W(y,x), 

(2) 

where da(y;mt,mt) is the leading-order partonic differ- 
ential cross section, q\ and q 2 are the momentum frac- 
tions of the colliding partons (assumed to be massless) 
within the incident p and p, and the sum runs over all 
possible combinations of initial-state parton flavors, jet- 
to-parton assignments, and all W — ► iv neutrino solu- 
tions In the sum over jet-to-parton assignments in 
P s ig, each permutation of jets carries a weight wi, which 
is the normalized product of probabilities for tagging any 
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FIG. 1: Values of the measured mean A from MC pseudo ex- 
periments as a function of A ln , parameterized by straight lines 
for (a) e+jets and (b) /i+jets MC events. Dotted lines rep- 
resent complete equality between measured and input values. 
Results from pseudo experiments with same A ln but different 
ra snm correspond to the extra points for fixed A m (see text). 



FIG. 2: Fitted contours of equal probability for the two- 
dimensional likelihoods as a function of mt and rxij for (a) 
e+jets and (b) ^i+jets data. The boxes, representing the bins 
in the two-dimensional histograms of the likelihoods, have ar- 
eas proportional to the bin contents, set equal to the value of 
the likelihood evaluated at the bin center. 



jet under a given parton flavor hypothesis The F(qi) 
include the probability densities for finding a parton of 
given flavor and longitudinal momentum fraction in the 
p or p assuming the CTEQ6L1 [Io| parton distribution 
functions (PDF), as well as the probability densities for 
the transverse components of the qi obtained from the 
LO event generator pythia llj . The normalization term 



tt 

norm 



The overall detection efficiency for tt depends on the 
values of both mt and mj. This is taken into account 
through the normalization by the observed cross section 
<orm = I A(x)P sig dx = a tt (mt,mt) (A(m t , m T )), where 
a u (m t ,mj) is the total cross section calculated by inte- 
grating the partonic cross section cr^| [13] , corresponding 
to the specific ME used in the analysis, over initial and fi- 
nal parton distributions and summing over initial parton 
flavors. (A(mt,rnj)) is the mean acceptance determined 
from the generated it events. The expressions for -Fkkg 
are similar, except that the probability does not depend 
on mt or mj. 

Samples of tt MC events with different values of and 
mj are required to simulate tt production and decay in 
order to calibrate the results of the analysis. These events 
are generated with a version of the PYTHIA generator 
modified to provide independent values of m t and mj. 
The specific values chosen for (mt-, mj) form a square 
grid spaced at 5 GeV intervals between (165,165) and 
(180,180), excluding the two extreme points at (165,180) 
and (180,165). The MC events for equal values of m t and 
mj are generated with the default version of pythia. 

Approximations made in formulating the likelihood 
can bias the final result. This issue is examined by com- 
paring the measured and input values of A in pseudo ex- 
periments composed of MC tt and VF+jets events. The 
calibration is shown in Fig. [T] in terms of the measured 
mean A as a function of its input value (A ln ), separately 
for the e+jets and /x+jets MC samples, for all MC sam- 
ples generated at the input reference points on the (mt, 



mj) grid. There are 2, 3, 4, 3, and 2 different (mt, mj) 
points with a common A ln of —10, —5, 0, +5, and +10 
GeV, respectively. The dispersions in the measured val- 
ues of mean A for different (m t , mj) points, but with 
same values of A in , are consistent with expected statis- 
tical fluctuations, as can be observed in Fig. [T] The fit 
X 2 /d.o.f. for the points in Figs. HJa) and QJb) are 1.8 
and 0.84, respectively. The parameterizations shown in 
Fig. Q] are used to calibrate L(x; A) for the selected data 
sample. 

We define the pull as (A - (A))/V(A) where A is the 
measured mass difference for a given pseudo experiment, 
(A) is the mean measured mass difference for all pseudo 
experiments, and tr(A) is the uncertainty of the mea- 
sured mass difference for the given pseudo experiment. 
The mean widths of the pull distributions for all sam- 
ples used in Fig. Q] arc 1.2 and 1.1 for e+jets and /z+jets, 
respectively. The deviations of these widths from 1 arc 
used to correct the measured uncertainties in data. 

Fitted two-dimensional Gaussian contours of equal 
probability (in terms of the standard deviation sd) for 
L(x; mt,mj) are shown for the electron and muon data 
samples in Figs. [SJa) and[2jb), respectively. The corre- 
sponding L(x; A) for both channels arc given in Figs.[3la) 
and [3jb). The two sets of data are consistent within 
their respective uncertainties, and the small correlations 
^e+jots = _ .05, ^+jct s = _q.01) extracted from the 
fits in Fig. [2] between m t and mj are not statistically sig- 
nificant, nor are the shifts in the projections shown in 
Fig. El 

Results from the two channels are combined through 
a weighted average of the separate electron and muon 
values. This has the advantage of using their respective 
pulls to adjust the uncertainties of each measurement 
before combining the two results. Using this averaging 
process, we quote the final combined means and their 
statistical uncertainties as A = 3.8 ± 3.4(stat.) GeV and 
m snm = 170.9 ± 1.5(stat.) GeV. The latter is consistent 
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A [GeV] A [GeV] 



FIG. 3: Projections of the likelihoods onto the A axis for (a) 
e+jets and (b) /i+jets data. 



with the previous measurement of Ref. [5| (see also Ref. 

The systematic uncertainties arc summarized in Ta- 
ble |H The first cate gory, Physics modeling, comprises the 
uncertainties in MC modeling of it and VT+jets events. 
The second category, Detector modeling, addresses uncer- 
tainties in the calibration of jet energy and simulation of 
detector response. The last category, Method, addresses 
uncertainties in the calibration and possible systematic 
effects due to assumptions made in the analysis. Ex- 
cept for two, all systematic uncertainties are identical 
to those described previously [f| . Many of these uncer- 
tainties (e.g., uncertainties in JES, PDF, jet resolution, 
multijet contamination) are expected to partially cancel 
in the measurement of the mass difference, but are of- 
ten dominated by the statistics of the samples used to 
evaluate them. The two new contributions address the 
possibilities of (i) reconstructing leptons with the wrong 
charge, and (ii) uncertainties from modeling differences 
in the response of the calorimeter to b and b jets [l4j . 
which can affect the measurement of the mass difference. 
These were evaluated for (i) by estimating the effect of 
an increase in charge misidentification in MC simulations 
that would match that found in data (~1% for both e 
and /z). For (ii), studies were performed on MC samples 
and on data seeking any difference in detector response 
to b and b quarks beyond expectations from interactions 
of their decay products, which are accommodated in the 
MC simulations. The observed differences were limited 
by the statistics of both samples. The total systematic 
uncertainty is 1.2 GeV. Combining the systematic and 
statistical uncertainties of the measurement in quadra- 
ture yields A = 3.8 ± 3.7 GeV, a value consistent with 
CPT invariance. 

In summary, we have measured the t and t mass dif- 
ference in ~1 fb _1 of data in £±jets it events and find 
the mass difference to be m t — mj = 3.8 ± 3.7 GeV, cor- 
responding to a relative mass difference of A/m sum = 
(2.2 ± 2.2)%. This is the first direct measurement of a 
mass difference between a quark and its antiquark part- 
ner. 



TABLE I: Summary of systematic uncertainties on A. 



Source Uncertainty (GeV) 

Physics modeling 

Signal ±0.85 

PDF uncertainty ±0.26 

Background modeling ±0.03 

Heavy flavor scale factor ±0.07 

b fragmentation ±0.12 

Detector modeling: 

b/light response ratio ±0.04 

Jet identification ±0.16 

Jet resolution ±0.39 

Trigger ±0.09 

Overall jet energy scale ±0.08 

Residual jet energy scale ±0.07 

Muon resolution ±0.09 

Wrong charge leptons ±0.07 

Asymmetry in bb response ±0.42 

Method: 

MC calibration ±0.25 

6-tagging efficiency ±0.25 

Multijet contamination ±0.40 

Signal fraction ±0.10 

Total (in quadrature) ±1.22 
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